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We analyse two issues that arise in the context of (matrix) string theories in plane wave back- 
grounds, namely (1) the use of Brinkmann- versus Rosen-variables in the quantum theory for 
general plane waves (which we settle conclusively in favour of Brinkmann variables), and (2) the 
regularisation of the quantum dynamics for a certain class of singular plane waves (discussing 
the benefits and limitations of regularisations of the plane- wave metric itself). 



1 Introduction 



As exact and potentially exactly solvable string backgrounds, plane waves have been extensively 
studied in string theory (see e.g. [1] for a review of the early literature). As such they also provide 
an ideal setting for studying the issue and fate of time-dependent backgrounds and space-time 
singularities in string theory [2, 3]. In the usual Brinkmann coordinates, a plane wave metric 
has the form 

ds 2 = 2dz+dz- + A ab (z + )z a z b (dz+) 2 + 5 ab dz a dz b , 

and this metric is singular iff the profile A a b(z + ) (= minus the matrix of frequency-squares for 
transverse string modes in the lightcone gauge z + ~ t) diverges for some (finite) value of z + } 
Among all singular plane waves, from the point of view of studying string theory in singular 
backgrounds, those with a profile of the form A ab (z + ) ~ {z + )~ 2 are singled out by at least two 
facts: 

1. Precisely these backgrounds were shown in [5, 6, 7] to arise universally as the Penrose 
limits [8, 9, 10] of a large general class of singularities (of power-law type [11]) which 
encompasses all the standard string theory black hole and black brane singularities as 
well as those of all (non-oscillatory/chaotic) standard cosmological models. Thus these 
special singular plane waves are not just toy-models of space-time singularities but honest 
lst-order approximations (in the spirit of the Penrose-Fermi expansion [12] around the 
Penrose limit metric) to realistic types of singularities. 

2. These singular plane waves possess precisely those symmetries (isometries and homoth- 
etics), in particular the boost/scale symmetry z^ — > X^z^, that are required [13] to 
implement the usual (flat space) Seiberg-Sen(-CSV) [14, 15, 16, 17] DLCQ derivation of 
matrix string theory [18] for these curved backgrounds, so that in principle one has avail- 
able a non-perturbative definition of string theory in these backgrounds (the plane wave 
matrix big bang models of [13]). 

The above-mentioned boost/scale symmetry renders the metric homogeneous away from the 
locus of the singularity at z + — (sec e.g. [10, 19, 20]). Thus these metrics can be considered 
as singular counterparts of the much studied (symmetric) plane waves with constant profile A ab 
(which are homogeneous because of the invariance under lightcone time translations), and this 
scale-invariancc is reflected in (and impacts upon) the quantisation of strings and particles in 
such a background [19, 20, 21, 22]. In the lightcone gauge z + = t, the string mode equations are 
simply harmonic oscillator equations with time-dependent frequencies 0J n {t) 2 = uj{t) 2 + n 2 , with 
uj(t) 2 ~ t~ 2 arising from the plane wave profile. In particular, perturbative string theory in these 
singular homogeneous plane waves (with weak string coupling at the singularity) was analysed 
in detail in [19], with the conclusion that, although there is divergent particle production on the 
world-sheet as one approaches the singularity at t = (cf. also [2] and [23] ) , there is a possibility 

1 What is meant by a singularity of the metric in the present context (where all scalar curvature invariants of 
the metric are zero) is divergent tidal forces, manifested by a divergence of components of the Riemann tensor 
in a parallel propagated frame (p.p.-singularities in the terminology of [4]). Thus plane waves with a singular 
profile have null p.p.-singularities. 
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that an analytic continuation of the string wave-functions could lead to non-singular evolution 
through the singularity. 

Curiously exactly the same (decoupled, Abelian) string mode equations arise from the analysis 
of strong string coupling singularities in the context of the (non- Abelian) matrix string models 
[13] for singular plane waves, in which the strong string coupling implies weak (time-dependent) 
Yang-Mills gauge coupling. This comes about because, as shown in [24], near the singularity the 
usual quartic interaction term for the non- Abelian matrix string coordinates becomes irrelevant. 2 
Thus the near-singularity behaviour of the matrix string coordinates is governed only by the 
decoupled (non-self- interacting) kinetic and harmonic oscillator (mass) terms ~ w n (i) 2 , as above, 
and one formally encounters the same problems and issues as in the perturbative string theory 
model of [19], now however at strong string coupling. 

The (thus necessary) issue of regularisation of the particle and string dynamics in these singular 
backgrounds was taken up again in [25, 26]. 3 In particular, in [26] a specific proposal was made 
to regularise the dynamics, namely by an explicit single-scale regularisation of the oscillator 
frequencies (equivalently, of the plane wave metric) of the form 

in which the scale invariance of the original background is as much as possible respected. In 
[26, 28] some minimal conditions were proposed for such a regularisation to lead to a well-defined 
quantum evolution in the limit e — > that the regulator is removed. One of the purposes of this 
article is to analyse and scrutinise this proposal in some more detail. 

Before embarking on this, however, there is another issue that needs to be addressed and settled 
once and for all which is a minor annoyance but which occasionally rears its ugly head. This 
is the question whether one should analyse or quantise a plane wave system in the Brinkmann 
coordinates used above, or whether one can (either alternatively or even profitably) study the 
system in variables based on Rosen coordinates, in which a plane wave metric takes the form 

ds 2 = 2dudv + gij{u)dx l dxi _ 

In particular, in the present context, metrics of this form appear to provide attractive models 
of null cosmologies [19]. However, as is well known (and as we will recall at the beginning 
of section 2), these Rosen coordinate systems are fraught with ambiguities, almost invariably 
exhibiting spurious coordinate singularities, while the Brinkmann coordinate system encodes 
directly geometrically invariant and physical information about the space-time (coordinates z a 
are measures of proper distance, and the profile A a b assembles the parallel orthonormal frame 
components of the Riemann tensor). 

Since it has, in spite of this, occasionally been suggested (see e.g. [29, 30]) that the Rosen system 
may provide a better, because seemingly less singular, description of the quantum dynamics, 
in section 2 we first quickly go through the (essentially undergraduate) exercise of constructing 

2 This is a stronger statement than just the fact that the Yang-Mills coupling constant goes to zero, its validity 
relying [24] on the precise time-dependence of the dilaton arising from the requirement that the metric-dilaton 
configuration be a solution to the string background equations of motion [13]. 

3 See also the review [27] for a discussion of the state of affairs. 
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and solving the quantum theory of Brinkmann and Rosen particles and the (unitary) relation 
between them. We then revisit the Brinkmann vs Rosen issue in this quantum context, stress 
that this is not a question of Rosen vs Brinkmann quantum states but rather a matter of which 
operators one attributes physical significance to, and show that the plane wave space-time origin 
of these models uniquely leads to the Brinkmann picture. This conclusion is then significantly 
strengthened and sharpened in section 2.4, where we show that some rather remarkable and 
striking recent uniqueness results [31, 32, 33, 34, 35] regarding the quantisation of scalar fields 
with time-dependent mass terms imply that for strings in arbitrary plane wave backgrounds (1) 
in the Brinkmann paramctrisation there is a unique Fock quantisation respecting the translation 
invariance along the worldsheet S 1 (namely the standard Fock space of a free massless scalar 
field) such that the dynamics is unitarily implemented; and (2) that there is no ^-invariant Fock 
quantisation of the Rosen system with unitary implementable time evolution (in the singular 
case, with a mass term ~ t~ 2 , these results hold for any t ^ 0). 

The upshot of this discussion is that the quantisation of (matrix) strings in plane wave back- 
grounds needs to be performed in Brinkmann variables (something that we had already advo- 
cated in [13] on the related grounds that it is only in these variables that the scalars have a 
canonical kinetic term) , and that in spite of the time-dependence of the mass-term / Hamilto- 
nian (which usually leads to quantisation ambiguities and phenomena like particle production 
etc.) there is a unique Fock quantisation with unitary time-evolution. 

Having disposed of the Rosen representation, we then (re-)turn to the issue of regularisation of 
the dynamics in singular scale- invariant plane waves and the proposal of [26] . We first specialise 
the results of section 2 regarding expectation values to the singular plane waves to illustrate 
how the plane wave singularity manifests itself in the quantum theory. We then exhibit a family 
of single-scale regularisations that satisfy the necessary conditions of [26, 28], essentially by 
a method that amounts to directly regularising the Rosen metric rather than the Brinkmann 
wave profile / frequency. We analyse the extent to which such regularisations can be considered 
to give a non-singular evolution through t = in the limit that the regularisation parameter 
e — > 0, with the conclusion that no regularisation of the single-scale form can provide this. More 
complicated regularisations that achieve this do exist, however, and we comment on this as well 
as on issues related to the regularisation of the dilaton. We end with a discussion of the results, 
in particular their implications for matrix string models of plane wave null singularities, and 
their interpretation in terms of the Penrose limit origin of these space-times. 

2 Quantum Brinkmann vs Rosen Plane Wave Dynamics 

In this section we will analyse the relation between the quantum dyanmics of particles and strings 
in the Brinkmann and Rosen parametrisations. To set the stage and motivate the discussion of 
this section, we begin by recalling the well-known ambiguities and spurious singularities afflicting 
the Rosen description of the classical particle dynamics. 
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2.1 Classical Dynamics and Rosen Singularities 



In Brinkmann (resp. Rosen) coordinates, plane wave metrics have the form 

ds 2 = 2dz+dz- + A ab (z + )z a z b (dz+) 2 + S ab dz a dz b 

. . (2.1) 
= 2dudv + gij(u)dx l dx J 

(u = z + , the remaining coordinate transformation between the unique Brinkman form of the 
metric and the non-unique Rosen form is given in (A. 3)). For present purposes nothing is lost by 
considering the isotropic (or transverse 1-dimcnsional) case gij (u) = e(t) 2 5ij and correspondingly 
A a b(u) = A(u)S ab . 4 The Brinkmann and Rosen Lagrangians for a relativistic particle in the 
lightcone gauge z + = u = t are 

L bc (z,z;t) = \{z 2 -u 2 {t)z 2 ) L rc (x;t) = \e(tfx 2 , (2.2) 

the general relation (A. 4) between the Rosen and Brinkmann parameters and A ab reducing 
to the harmonic oscillator equation 

A(t) = ^ = -cu 2 (t) or e{t)+io 2 {t)e{t) = . (2.3) 
e(t) 

Up to a total derivative term, L bc and L rc are related by the coordinate transformation z = e(t)x, 
the former being exactly the Lagrangian of a time-dependent harmonic oscillator, the latter 
that of a "free" particle with time-dependent mass. The non-uniqueness of Rosen coordinates, 
x = e(t)~ 1 z, and of the Rosen Lagrangian, is reflected in the arbitrary choice of solution e(t) of 
the oscillator equations of motion (2.3). The corresponding Hamiltonians, 

H bc = \{p 2 z + io(t) 2 z 2 ) H rc = \e(t)- 2 p 2 x , (2.4) 

are then related by a linear time-dependent canonical transformation, the only (barely) non- 
trivial feature being that as a consequence of this time-dependence the Hamiltonians differ not 
only by the canonical transformation of the canonical variables but also by the time-derivative 
of the corresponding generating function (A. 8). 

Brinkmann coordinates provide a global coordinate chart for plane wave space-times while Rosen 
coordinates typically (and almost invariably, e.g. for backgrounds satisfying the weak energy 
condition, see [36] and Appendix A.l) exhibit spurious coordinate singularities. As a prototypical 
example, consider the case of constant A ab i.e. constant frequency oj. In Rosen coordinates one 
attempts to use the classical trajectories of the Brinkmann system as Cartesian coordinate lines 
(as reflected in the fact that x k = const, is always a solution of the Rosen equations of motion), 
which in the present case obviously doesn't bode well for these coordinates as the oscillator 
reaches its turning point. Indeed, with the choice e(t) = sinwi, say, the Rosen metric and 
Hamiltonian have the form 

ds 2 = 2dudv + sin 2 Lot dx 2 H rc = i(sin 2 ^t)~ x p 2 x (2.5) 

4 In particular, such metrics can always be elevated to full-fledged string theory background by a suitable 
choice of null dilaton. 
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which is evidently not able to describe the perfectly regular evolution of the system (or dynamics 
of strings in e.g. the maximally supersymmetric IIB background [37, 38, 39, 40, 41]) beyond 
t = tt/uj. 

The situation may at first appear to be less clear-cut in the case of genuinely singular plane waves 
in which both Rosen and Brinkmann metrics exhibit a singularity at, say, u — 0. As mentioned 
in the introduction, Rosen metrics are occasionally regarded as null analogues of cosmological 
metrics. For instance, isotropic metrics with gij{u) = g{u)5ij look like null counterparts of 
standard spatially flat {k = 0) FRW metrics, and metrics with a power-law behaviour gij(u) = 
u 2ai 6ij are null analogues of Kasner metrics. The latter, incidentally, are precisely the Rosen 
forms of the scale-invariant plane waves for a certain range of frequencies, 

ffy(u) = u 2a *% => A lJ (z+)=a l (a l -l)S lJ (z+y 2 . (2.6) 

However, interpreting Rosen plane waves in this way requires some care. For example, consider 
the simple (both isotropic and power-law) metric with g(u) = u 4 . This form of the metric 
suggests that it describes a null version of a Big Bang at u = 0. However, the very same metric 
can be written in a different coordinate system as a Rosen plane wave with g(u) = u~ 2 , with 
the same u, suggesting that what happens at u — is a Big Rip, not a Big Bang. Clearly 
then, the connection beween the form of the metric and the physics and geometry it actually 
describes is at best somewhat hidden in Rosen coordinates. It is however completely manifest in 
Brinkmann coordinates (in which e.g. both the above models give rise to a scale-invariant plane 
wave with uj 2 {t) = -It- 2 ), since Brinkmann coordinates are Fermi coordinates along the family 
of null geodesies with z a = and z~ = const. [12]. In particular, the transverse coordinates z a 
are Riemann normal coordinates and thus a direct measure of physical proper distance in the 
space-time, and the wave profile A a b(z + ) = —R +a+ i, gives the (only non-vanishing) components 
of the Riemann tensor in a parallel propagated frame along this null geodesies, and thus directly 
encodes the only non-trivial geometric features of a plane wave in a generally-covariant way. We 
will come back to this issue in the context of the corresponding quantum systems in section 2.3. 

2.2 Rosen and Brinkmann Quantum States 

For these simple classical systems, the quantum dynamics is evidently determined completely 
by the classical dynamics, i.e. can be expressed in terms of solutions to the classical equations 
of motion. At first this seems to favour the Heisenberg picture, but it turns out that, in 
order to disentangle the time-dependent canonical transformation from the time-evolution, the 
Schrodinger picture is marginally more convenient and provides some more insights into the 
properties of the Rosen quantum system. Moreover, the proposal of [25], which is the subject 
of section 3, is phrased in terms of the (WKB-exact) WKB wave function of the Brinkmann 
system. For both these reasons, we will work in the Schrodinger picture. 

As generically both the Brinkmann and the Rosen Hamiltonians are time-dependent, we will 
then primarily be interested in the solutions to the time-dependent Schrodinger equations 

idtM*) = H bc ipt(z) id t i>' t {x) = H rc ^ t (x) (2.7) 
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and in the time-dependence of expectation values of suitably chosen operators. Since the canoni- 
cal transformation between the Brinkmann and Rosen systems is linear (and we are dealing with 
a finite-dimensional quantum system) , on general grounds there exists a unitary transformation 

U : Urc -> H bc : Uijj 1 = V (2.8) 

from the Rosen to the Brinkmann space of states implementing this canonical transformation 
at the quantum level and mapping solutions of (2.7) into each other. 5 Its explicit form will not 
be needed in the following but is given in Appendix A. 2 for completeness' sake. 

Turning now first to the Rosen system, the Schrodinger equation is 

idti/Ux) = Hrcttix) = ^LfftMx) (2.9) 

In terms of the Rosen conformal-time coordinate T(t) (see Appendix A.l), (2.9) reduces to the 
Schrodinger equation of free particle in flat space, 

dT = dt/e{tf id T ip' T (x) = \y\^ T (x) . (2.10) 

Thus a convenient and simple basis of solutions to (2.9) is provided by the Rosen momentum 

^ k{x)=e -(^/2)T(t) + t kx (2U) 

(these are evidently eigenstates of p x with eigenvalue k, delta-function normalised in the usual 
way). Note that the fact that the transformation t — > T(t) maps the Rosen model evolution to 
that of a free particle in flat space already shows that in general such a description cannot be 
valid globally (and this reflects the Rosen coordinate singularities). 

The time-dependent Schrodinger equation of the Brinkmann system 

idtM*) = HbcMx) = Wl + "(tfz 2 )Mz) , (2-12) 

on the other hand, is simply that of a time-dependent harmonic oscillator. Since we have a free 
(quadratic) theory, the path integral can be explicitly calculated, 



Here z c (t) now denotes the classical harmonic oscillator solution with the given boundary con- 
dition, z c (ti) — Zi and z c (tf) — Zf, S[z c ] is the classical action, and the fluctuation determinants 
are to be calculated for zero (Dirichlet) boundary conditions. 

Both the 1-loop fluctuation determinant and the classical action can be compactly expressed in 
terms of the Gelfand-Yaglom (GY) function (see e.g. [42] and references therein), the solution 
of the classical equations of motion characterised by 

{d 2 t +u(t) 2 )F u (t) =0 F ti (U) = F tt (U) = l (2.14) 



5 The notation throughout will be such that wave functions in Brinkmann coordinates (Brinkmann states) will 
be denoted by ■tfj(z), perhaps with other labels indicating individual states or a basis of states, while Rosen states 
will be denoted by xp'(x). Wc will not keep track of factors of 2tt in the normalisation of the wave functions. 
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(and its dual G ti (t) — —d ti F ti (t), a linearly independent solution characterised by G ti {ti) = 1, 
Gti{ti) — 0, with Wronskian W(F, —G) = +1). In terms of any linearly independent pair e\^{t) 
of the classical equations of motion, the GY solution can be constructed as 

F t .(t) = ^Ml^M) . (2.15) 
W(ei,e 2 ) 

For the classical action one then has the expression 

S[z c } = ^^[F u (t)z 2 + G ti (t)zf-2zz i ] , (2.16) 
and for the ratio of fluctuation determinants one has the remarkably simple result 



1 Det[-9 t 2 ] 1 



j2mh(t f -U)\ Det[-d t 2 - u{tf] ^2nihF ti {t f ) 



(2.17) 



The above-mentioned properties of the GY function and its dual imply that it also governs the 
time-evolution of the classical system, since one has 

(z(t)\ = (G u {t) F u (t)\ (z{U)\ , . 

\m) \G u (t) Kit)) \z{u)) 1 • > 

The amplitude (2.13) therefore gives rise to the family of Brinkmann propagator (or WKB) 
states 

ip tiZ .(z) -i P iS[z c (z,Zi;t,ti)} ( 219 ) 

V F u (*) 

characterised by lim t _,. ti ipt, Zi (z) = 5{z — Zi). Thus the GY function compactly encodes much 
of the relevant information about the quantum system, and for this reason we will focus on the 
behaviour of the GY function (and expectation values of certain operators) in our analysis of 
singular systems in section 3. 

2.3 Brinkmann vs Rosen expectation values 

Turning now to expectation values of quantum operators, a particularly simple (and easy to 
work with) basis of states is provided by the Rosen momentum states ip' t k (2.11). Thus we 
consider general wave packets 

i// t {x) - J dk a(fc)$, fc (aO (2-20) 

(with the assumption that the a(k) are square integrable and differentiable) and their Brinkmann 
duals ip = Utp'. For position observables, Rosen and Brinkmann expectation values are con- 
nected in a simple way, U~ 1 zU — e(t)x (A. 12) implying relations like 

<^\z 2 \4,>=e(t) 2 <4>'\x 2 W > . (2.21) 

Moreover, since the Rosen momentum states are eigenstates of the momentum operator p x , it 
is straightforward to determine the expectation value of the Rosen Hamiltonian in such a wave 
packet, 

< ij'\HrcW >= J dk k 2 a{k)*a{k) . (2.22) 
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We note that if the wave packet ip' t x is to be a solution of the Schrodinger equation, the coef- 
ficients a(k) need to be time-independent. Thus we conclude that the time-dependence of the 
expectation value of the Rosen Hamiltonian is determined entirely by the prefactor e(i)~ 2 of 
the Hamiltonian. In particular, if e(t) has a zero somewhere, e(to) = 0, the expectation value 
of the Rosen Hamiltonian diverges as t — > to f° r an y solution to the Schrodinger equation. We 
will come back to this and related issues below. 

In order to calculate (and compare this with) the expectation value of Hb c , one needs to take 
into account the non-trivial relationship between the classical Hamiltonians provided by the 
generating function (A. 8). Relegating the details of the calculations to the appendix, the results 
for the expectation values of H rc , x 2 and their Brinkmann counterparts can be summarised as 

(ip'\x 2 \il>') = A + T(i)C + T(t) 2 M 
(ip\z 2 \ip) = e(t) 2 A + e(t) 2 M + e(i)e(t)C 
{il>'\HrcW) = \e{t)- 2 M 
(ip\H hc \Tp) = E(e)(t)A + E(e)(t)M + E(e, g)(i)C 



(2.23) 



where the (constant) coefficients A^0,B^O,C, defined in (A. 22), depend only on the modes 
a(k) of the wave packet, and the functions 

£(e) = \{<? + ^ 2 e 2 ) and E(e) = ^(e 2 + u?e 2 ) (2.24) 

can be readily intepreted as the classical energies of the solutions e(t) and e(t). We have also 
introduced E(e, e) = l/2(ee + w 2 ee) for uniformity of notation, and we have used some identities 
satisfied by e(t) and its dual linearly independent solution e(t) = e(t)T(t) to put the results for 
the Brinkmann expectation values into a form which makes it manifest that it is symmetric 
under the exchange e(t) o e(i). 6 

None of the results (2.23) should come as the slightest surprise since they just mirror the time- 
dependence of the corresponding general classical solution, a linear combination of e(i) and its 
dual e(t) = e(t)T(t). In particular, this gives the expected result for the expectation value of 
the Hamiltonian of the usual time- independent harmonic oscillator with constant frequency to. 
In this case one has e.g. e(t) = sinwi and e(t) = — uj^ 1 cosujt so that 



e(t) = sinwt E(e) = \u? , E(e) = \ , E(e,e) = 



(2.25) 



For the ground state wave function, with a(£) ~ exp(— £ 2 /2u), one finds (il)\Hb c \if>) = to/2. 

Note that the corresponding Rosen Hamiltonian, on the other hand, fails completely to exhibit 
this perfectly regular behaviour of the time-independent harmonic oscillator. As noted before, 
the expectation value of the Rosen Hamiltonian diverges for all states as t — >• to with e(to) =0, 
and for any solution of the harmonic oscillator equations there is an infinite number of such to. 



6 Recall that, unlike the Rosen system, the Brinkmann system depends only on the frerquency ui(i), not on a 
particular classical solution e(t), and this should be (and is of course) reflected in the expectation values etc. 
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Geometrically, this is a pure coordinate singularity of the Rosen metric, and the quantum system 
fails to detect that this is not an honest singularity. The problem is not so much the Rosen states 
but the operators whose expectation values one wants to evaluate. To highlight this, consider 
the statements 

<^|ff 6c |V>> = <^'|(#6c)rc|V''> < V\HrcW > = <M{H r c)bc\^> (2.26) 

that follow from (A. 16). 7 In the time-independent case, say, of a harmonic oscillator with 
constant real frequencies, the Brinkmann expectation value of the Brinkmann Hamiltonian is 
of course well-defined and non-singular in every Schrodinger state. Therefore, by the above 
identity, also the expectation value of the corresponding Rosen operator (Hbc)rc is non-singular 
and well-defined in every one of the corresponding Rosen states - it is a good Rosen operator. 
The Rosen Hamiltonian H rci on the other hand, is not: it has a divergent expectation value in 
any Rosen state, and its Brinkmann-picture counterpart (H rc )bc correspondingly has a divergent 
expectation value in every Brinkmann state. 

The fact that the Brinkmann versus Rosen issue is a question of operators and not of states is also 
reinforced by the space-time origin and interpretation of these models, which provide additional 
physical criteria allowing one to select the operators one attributes physical significance to, thus 
going beyond the mere quantum mechanical formalism. 

As an example, let us consider the Rosen expectation value < x 2 >. This gives the mean value 
of the Rosen coordinate-distance squared in some state. Without any further physical input, 
one can argue until one is blue in the face whether this is the physically relevant quantity to look 
at or whether one should rather look at < z 2 >. Indeed, in principle nothing stops one from 
constructing a non-relativistic quantum system with a Rosen Hamiltonian, and with the Rosen 
coordinates x l Cartesian flat space coordinates, and in this case indeed it would be < x 2 > that 
is the relevant quantity that can be related to observations and measurements of the system. 

However, recalling the origin of the Rosen and Brinkmann Hamiltonians as the lightcone Hamil- 
tonians of relativistic particles in the plane wave background, albeit in different coordinate 
systems, the situation presents itself in a rather different light. Indeed, in this case physical sys- 
tems detect and measure not coordinate distance but generally-covariant quantities like proper 
distance, and this can easily be accounted for, 

< (proper Rosen distance) 2 >=< e(t) 2 x 2 > . (2.27) 

But according to (2.21), this is just the Brinkmann expectation value of z 2 , 

< V'' | (proper Rosen distance) 2 ^' >=< ip\z 2 \ip > . (2.28) 

In this case, even if one starts in the Rosen system one is invariably led to the Brinkmann 
expectation value as the physically well-defined quantity to look at. 

Since the time-evolution of this physically meaningful notion of distance is governed by the 
Brinkmann Hamiltonian and not by the Rosen Hamiltonian (dz/dt ~ [Hb c ,z]), this shows that 

7 (Hb c )rc denotes the Brinkmann Hamiltonian expressed in terms of the Rosen operators x,p x etc. Note that 
this is not the same as the Rosen Hamiltonian H rc , cf. (A. 17) for the relation between the two. 



10 



it is also the Brinkmann and not the Rosen Hamiltonian that is a physically meaningful operator 
in the quantum theory of a relativistic particle. 

We will see below that this conclusion can be singificantly strengthened in that the Brinkmann 
dynamics is also singled out in the dynamics of strings and fields by the requirement of unitarity 
(and that the Rosen dynamics cannot be unitarily implemented on any field theory Fock space) . 

2.4 Unitarity and Uniqueness at the QFT/String Level 
(Brinkmann wins hands down) 

So far we have studied aspects of the quantum mechanics of point particles in plane wave 
backgrounds. When one considers perturbative string theory, the relevant dynamics of closed 
strings (in the lightcone gauge) is that of the transverse string modes on the cylinder S 1 x R, 
and it is well known that in Brinkmann variables this reduces to the standard free action of 
massive scalar fields with time-dependent mass matrix — A a b{t) (the frequency /profile of the 
Brinkmann plane wave). The same is true for non-perturbative matrix string theory: in [24] it 
was shown that the near-singularity t — > evolution of the plane wave matrix big bang models 
of [13] is dominated by the same time-dependent mass terms, while the quartic interaction and 
the accompanying non-Abelian nature of the dynamics are subleading. In the isotropic case we 
have focussed on one is thus in both cases confronted with the problem of the quantisation of a 
single scalar field on S 1 x M with time-dependent mass term m 2 (t) = cj 2 (t) (or with an infinite 
number of Fourier modes with masses m n (t) 2 = m(t) 2 + n 2 , n e Z). 

It is well known that in the quantisation of systems with an infinite number of degrees of freedom 
(for which is there is no analogue of the Stone - von Neumann theorem) one is confronted with 
various ambiguities, e.g. in the choice of field parametrisation or, for a fixed parametrisation, in 
the choice of (Fock) representation of the canonical commutation relations (CCRs). In the usual 
setting of a Poincare-invariant field theory (of a free scalar field, say), the requirement of Poincare 
invariance of the Fock vacuum selects a unique Fock representation of the CCRs on which the 
dynamics can be unitarily implemented. In general (in particular non-stationary) space-times 
there is no such unique choice, and this give rise to phenomena like particle production etc. 

When there are some spatial isometries, it is natural to require that the quantisation procedure 
respect these symmetries, but (a) even in spatially maximally symmetric (cosmological) situa- 
tions this is not enough to uniquely fix the representation, and (b) this says nothing about the 
dynamics and the possibility to unitarily implement the time-evolution. 

In a series of papers, motivated by cosmological considerations, Cortez et al. (see e.g. [31, 32, 33, 
34, 35]) investigated whether and to which extent the joint conditions of (a) invariance under 
the background symmetries and (b) unitary implementation of the dynamics can be used to 
select and specify a unique Fock quantisation. In particular, they investigated the case of scalar 
fields on S d x M, d < 3, with a time-dependent mass (precisely the situation we are interested 
in here for d = 1), with some rather remarkable and striking uniqueness results. 

To describe these, consider a scalar field (f> satisfying the equation 

4> - Ac/) + s(t)<t> = (2.29) 
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where A is the Laplacian on S d and s(t) is the mass-squared function satisfying some mild reg- 
ularity condition (C 1 on the time interval / under consideration being sufficient). In particular, 
no condition on the sign of s(t) is (or needs to be) imposed. 8 

SO (d + l)-invariance is now imposed as a condition on the complex structure defining the 
creation and annihilation operators (a„,d*), thus on the Fock representation of the CCRs (this 
still leaves a lot of possibilities). The requirement of unitary implementability of the dynamics 
is then tantamount to the condition that the Bogoliubov transformation 



describing the temporal evolution U — > t of some suitably defined instantaneous creation and 
annihilation operators leads to finite "particle production", ^2 n \(3 n \ 2 < 0. Then the first main 
result is the following: 

1. For every s(t) (subject to the above-mentioned regularity assumptions), up to unitary 
equivalence there exists a unique SO(d + l)-invariant Fock representation of the CCRs 
(complex structure) with respect to which the dynamics is unitarily implemented. This 
representation is independent of s(t) and, in particular, coincides with the standard Fock 
representation of a free massless field in Minkowski space. 

Note that spatial compactness is crucial for the validity of the 2nd statement, IR problems 
preventing an analogous statement from being true e.g. for a free scalar field in Minkowski 
space, where the Fock represenations (and corresponding representations of the Poincare group) 
are of course not unitarily equivalent for different masses. 

In the above, we started with a fixed parametrisation of the fields, the scalar field <p an d its 
canonically conjugate (a spatial density on S d ). The second main result of Cortez et al. [35] 
concerns the ambiguity resulting from time-dependent linear canonical transformations (field 
redefinitions) of the form 



where h is the determinant of the spatial metric, F{t) is different from zero everywhere and 
both F(t) and G(t) are differentiable. Such transformations are commonly performed e.g. in 
the context of cosmological perturbation theory where they are used to put the action for the 
perturbations into a suitable (canonical) form. This transformation has an obvious and simple 
effect on the above Bogoliubov transformation implementing the time evolution, and one can 
analyse the conditions arising from the square-summability of the new \f3 n \. It turns out that 
the joint requirements of symmetry and unitarity also determine an essentially unique field 
parametrisation: 

8 A caveat in e.g. [33, 35] regarding such a condition for the quantisation of the zero mode sector is unnecessary; 
the quantum theory of a particle with potential s(t)z 2 , say, is well-defined for any sufficiently regular s(t), 
regardless of its sign; one should simply not try to define it via naive analytic continuation from the usual 
Schrodinger theory of a harmonic oscillator with u) 2 > as that would lead to a non-hermitian Hamiltonian. 




(2.30) 



<P = F(t)<t> 



P v = P*/F(t) + G(t)Vh<t> 



(2.31) 
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2. The dynamics cannot be unitarily implemented in an SO(d + l)-invariant Fock represen- 
tation unless F(t) and (for d ^ 1) G(t) are constant. 

When d = 1, shifts P^ — > P^ + G(t)<fi of the momentum are allowed, all leading to the same 
unique (free massless) representation of the CCRs singled out in the first result. We will see 
below that, for the point transformation that arises when going from Brinkmann to Rosen 
coordinates, F(t) constant already all by itself implies G(t) — so that this is, in any case, not 
an issue. 

Intuitively, the fact that there is a unique field parametrisaton that allows a unitary imple- 
mentation of the dynamics can be attributed to the fact that any non-trivial time-dependent 
reparametrisation of the field will give rise to damping/friction terms in the equations of motion 
which are obstructions to unitarity [35]. 

These results have several implications for the study of (matrix) string theory in plane wave 
backgrounds, in particular as regards the issue of whether to use a Brinkmann or a Rosen 
parametrisation of the fields, in which the scalar sector of the action schematically takes the 
form 

S bc = -3 / d2(J (v aP d a Z a d f3 Z b - A ah (t)Z a Z b + ...) (2.32) 

or 

Src = -\J d 2 a {^gaffydaJCdpX 1 + . . .) , (2.33) 

the (+ . . .)-term including possible interaction terms and/or couplings to background fields. 
Exactly as for the point particle action (Appendix A.l), up to a total derivative these two 
actions are related by the time-dependent scaling X 1 — e l a (t)Z a with gij(t)e l a (t)el(t) = 5 a b 
(this equivalence actually extends to the full non-Abelian matrix string or 3-algebra multiple 
M2-brane action [43]). In the isotropic case g%j{t) = e(t) 2 Sij, this corresponds to the canonical 
transformation (A. 7) 

X = Z/e(t) , P x = e{t)P z - e(t)Z , (2.34) 

which is precisely of the form (2.31) with F(t) — l/e(£) and G(t) = — e(t) (in particular, in this 
case one has F(t) = const. G(t) = 0). Note that the Brinkmann field Z has a canonical 
kinetic term and therefore satisfies the canonical equations of motion (2.29), while the Rosen 
field X plays the role of the reparametrised field ip. 

The implications of the results of Cortez et al. are now clear. We consider frequencies / mass 
terms s(t) = oo(t) 2 = —A(t) and a time interval / such that the assumptions about s(t) are 
satisfied (and will comment on the implications for more singular situations in section 3). 

1. In the Brinkmann parametrisation there is a unique Fock quantisation respecting the 
translation invariance along the worldsheet S 1 (namely the standard Fock space of a free 
massless scalar field) such that the dynamics is unitarily implemented. 

2. There is no ^-invariant Fock quantisation of the Rosen system with unitary time evolu- 
tion. 
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A priori these results do not exclude the possibility that, starting from scratch in Rosen coor- 
dinates, one could find an exotic Fock representation of the CCRs and implement unitarily the 
temporal evolution. We know for sure, though, that the vacuum of such a representation would 
not be invariant under S 1 rotations. The other alternative would be a non-Fock representation. 
Neither option seems palatable. We believe that this should settle once and for all the question 
whether one should use Brinkman or Rosen fields, and we will focus entirely on the Brinkmann 
systems in the following. 



3 Regularisation of the Dynamics for Singular Plane Waves 

We will now consider the quantum models arising from Brinkmann plane wave metrics (A. 2) 
with profile A a b(z + ) <~ (z + )~ 2 , as these were shown in [5, 6, 7] to arise generically as the Penrose 
limits of metrics with singularities, and to which the derivation of [17] of a matrix string model 
(the matrix big bang) was extended in [13]. The frequencies that typically arise from the Penrose 
limit procedure have the form a(l — a)(z + )~ 2 for some a £ R and continuing to focus on the 
isotropic case we will thus specifically study the models with frequency / mass terms 

o , \ ci(l — a) , 
w (*) = ^T 2 - ■ 

The problem is evidently invariant under a — > 1 — a and we can therefore, without loss of 
generality, restrict to the range a > 1/2. 9 We will mainly be interested in the range a > 1 
for which ui 2 (t) is negative, since it is this tachyonic behaviour that characterises the Penrose 
limits of string theory backgrounds with strong string coupling singularities [13] (for which a 
non-perturbative description such as matrix string theory thus also becomes mandatory), but 
many of our considerations are also valid in the weakly coupled range 1/2 < a < 1 which was 
the focus of interest e.g. in [19]. 

We recall and reiterate here that it was argued in [24] , on the basis of numerical and analytical 
investigations of various toy-models, that both the classical and the quantum evolution of the 
matrix (string) system near strong string coupling (and hence weak gauge coupling) singularities 
are generically driven entirely by the divergent tachyonic mass terms. In particular, the quartic 
interaction potential, which comes with a time-dependent coefficient ~ t 2q with q > 0, 

can be shown to be self-consistently small in quantum-mechanical perturbation theory in that 
regime and appears to (perhaps somewhat disappointingly) play no role there. Thus, in the 
following we can and will focus on the linear harmonic oscillator dynamics even when what we 
have in mind arc primarily applications to matrix string theory. 



3.1 Necessity of Regularisation 

The considerations of section 2.4 imply that there will be a unique Brinkmann quantisation of 
the system with a unitary time evolution as long as one stays away from t = (since (3.1) 

9 The value a = 1/2 would require a separate treatment, because one of the classical solutions has a logarithmic 
~ t 1 / 2 logt rather than a power-law behaviour, but is not particuarly interesting in other respects. Aspects of 
the a = 1/2 model (referred to as the Gowdy model in the cosmological context) are analysed in detail in [31]. 
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satisfies the required regularity assumptions for t > to > 0). 

On the other hand, it is self-evident that the quantum theory with frequency / mass terms will 
not miraculously regularise the classical singularity of the theory as t — > 0, since the evolution 
of wave packets and expectation values of this WKB-exact quantum system will track that of 
the singular classical system. Nevertheless, we will now explicitly display some of the resulting 
singularities, for reference purposes, and since we will attempt to regularise these expressions 
below. 

Starting from the fact that the two linearly independent solutions of the classical equations of 
motion 

e(t) + a(l - a)t~ 2 e(t) = (3.2) 
of the point particle system, defined for t > 0, can be chosen to be 

e(t) = t x - a , g(t) = t a /(l - 2a) (3.3) 

(with unit Wronskian), one finds from (2.15) that the Gelfand-Yaglom function is 

F * i (*) = T32^(^ 1 " a -*i"°* a ) ■ ( 3 - 4 ) 

A priori this is not defined for ti or t negative. Moreover, this expression evidently diverges 
for all t > as U — > for a > 1 and goes to zero for all t > as ti — > for 1/2 < a < 1, 
either behaviour indicating a breakdown of the quantum theory. Note that the appearance of 
an isolated zero of the GY function for some t = tf, F ti (tf) = 0, is not a fundamental problem. 
In view of (2.14), this just signals the existence of a zero mode of the Dirichlet problem which 
needs to (and can) be taken care of in any number of standard ways. What happens here is 
that the GY function goes to zero (or diverges) for all t > 0, and this is a genuinely singular 
behaviour. 

One can also see that unitarity of the time-evolution breaks down as t — > 0, signalled by the 
divergence of the Bogoliubov coefficient (3 n =o{t,U) m (2.30) as ti — > 0, leading a fortiori to a 
breakdown of unitarity of the full quantum field / string theory. This can be established either by 
transforming the GY time-evolution matrix (2.18) to an appropriate oscillator basis, or directly 
by a calculation of the Bogoliubov coefficient j3(t, U) between instantaneous eigenstates at times 
t and U. w 

Finally, for Brinkmann expectation values one finds from (2.23) that schematically 

(ip\z 2 \ip) ~ At 2 - 2a + Mt 2a + Ct 

(3 5) 

(ilj\Hbc\ip) ~ At~ 2a + mt 2a ~ 2 + ct~ 1 , 

which are also all singular as t — > (suppressed numerical factors hidden in the ~ signs making 
the coefficients of the potentially divergent terms zero only for a = 0, 1, i.e. uj 2 (t) = 0). 

"For a ^ 0, 1, /3(t,U) turns out to have the general form "(terms ~ t- a ,t 1 - a ) X (terms ~ t"~\*f)-(« <-> U)" . 
Thus as ti — > 0, 0(ti, (2) diverges for any a ^ 0, 1 (for a = 0, 1 one just has a free particle). 
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3.2 Single Scale Regularisation, Gelfand-Yaglom and Fine-Tuning 



We will now try to make (at least some of) the above expressions well denned by a regularisation 
of the frequency / plane wave profile (thus maintaining the (a'-exact) plane wave character of 
the background metric). Adopting at first the suggestion in [26] we will explore the possibility 
of using a single-scale regularisation of the form 

L0 2 {t) -> n 2 Jt) = e- 2 n 2 (t/e) = e- 2 n 2 (f]) (3.6) 

in which the scale invariance of the original background is as much as possible respected and 
the classical equation of motion becomes independent of e, 

^f(r 1 ) + n 2 (r ] )f(r ] )=0 . (3.7) 

In order to qualify as a regularisation, S1%(t) should be non-singular for all t, in particular also 
t = 0, as long as £ ^ 0, and should reduce to uj 2 (t) as e — » for t > 0. In terms of Q 2 {rj) this 
amounts to the conditions 

3 lim Q 2 (r]) = O 2 (0) and ft 2 (77) -S- ct(l - a)^ 2 for 77 -> 00 (3.8) 

7)— >0 

We will moreover assume that for e — > and i 7^ one reaches the same limiting frequency 
a(l — a)t~ 2 for t < as for t > (the case of different left/right frequencies having already 
been dismissed in [26]). In that case we also have the same asymptotic behaviour of ft 2 (r?) as 
77 — > —00. 

Some examples of rcgularisations satisfying these conditions that come to mind immediately are 
n 2 (t) = a(l - a)(t 2 + e 2 )- 1 & Q 2 (r?) = o(l - a)( V 2 + l)" 1 (3.9) 
and its generalisations 

tf( V ) = a(l-a) {r)2 l i)n+1 (3.10) 

(all of which we will unfortunately have to dismiss below for not accomplishing what they are 
supposed to do). Note that these examples have the property that fi 2 (— 77) = il 2 (77) , so that 
the regularised frequency preserves the Z 2 reflection symmetry t — )• — t of the unregularised 
frequency. This is a natural additional assumption, but one which will not play a significant 
role in the analysis below. 

We can then assume (see also [28]) that we have two linearly independent and e- independent 
solutions fk(v) °f (3-7) with the asymptotic behaviours 

(3.11) 

(Wronskian conservation implying that the ratio of the determinants of the two matrices is —1). 
Note that the generic behaviour of a solution as |?7| — > 00 is the dominant behaviour / ~ |7i| a 
(recall that we have chosen to consider the range a > 1/2). We will call a solution fine-tuned if 
it approaches the subleading solution | ^7 1 1 tt for both r\ — > +00 and 77 — > —00. The existence of 
such a solution is not guaranteed. If it exists, it has the following properties: 
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• If a fine-tuned solution exists, it is unique (any other linearly-independent solution has to 
exhibit the dominant asymptotic behaviour ~ |7y| a for a non-zero Wronskian at ±00). 

• If the regularised potential is even, then the fine-tuned solution is either even or odd (since 
otherwise one could construct another fine-tuned solution by taking its even or odd part). 

From this pair of linearly independent solutions we can constuct the Gelfand-Yaglom function 

Fl(t) = Wifuh)- 1 (fi(U/e)f 2 (t/e) - f 2 (U/e)f 1 (t/e)) (3.12) 

and we want to investigate its behaviour as e — > 0. By construction this limit is well defined if 
ti,t are both positive, giving rise to (3.4). The main result is that the limit is also well defined 
for ti < 0, t > iff a fine-tuned solution exists, 

U<0,t>0: 31im F t e .(t) 3 fine-tuned solution (3.13) 

Proof: in evaluating this limit using the asymptotic behaviour, and keeping in mind that a 
constant 77- Wronskian W v (fi,f 2 ) implies that the relevant t- Wronskian, t = er], is <~ e _1 , one 
finds one potentially divergent contribution <~ e 1 ^ 2 " whose coefficient is a+7- — a_7 + . But 

a+7_-a_7+=0 <S> 3A, C : { Aa ++ C ^+^° ^ Afi+Cfc fine-tuned (3.14) 

I Aa_ + Cj- = 

This is a reformulation of a statement in [28, 26] (where the GY function is called the compression 
factor). With this condition in place, one finds that the regularised GY function is finite as e — > 0. 
Assuming without loss of generality that in (3.11) it is f 2 that is the fine-tuned solution (so that 
7 ± = and a + a_ 7^ 0), one finds 

limF ti (*) = —!— (q\t\ a \t i \ 1 - a + q-%\ a \t\ 1 - a ) , (3.15) 

e^O 1 — la 

where q = ol+/ol—. Whenever the fine-tuned solution f 2 can be chosen to be of definite parity 
(e.g. when the regularised profile is even), one has |g| = 1 (q = ±1 for f 2 odd/even). Note 
that this is quite different from a naive analytic continuation of (3.4) through t = 0, which (for 
non-integer a) would give a relative imaginary phase between the two terms, and from the result 
one obtains [19] from the analytic continuation of the Besscl functions appearing as solutions to 
the string mode equations in this background. 

This raises the questions how restrictive the condition is that the regularised frequency admits a 
fine-tuned solution and how to characterise the regularised profiles that admit such a fine-tuned 
solution. We will now provide some partial answers to these questions. 

We consider the range a > 1 for which uj 2 (t) = a(l — a)t~ 2 < 0. First of all, we observe that 
a fine-tuned solution, behaving as |r?| 1_a for \rj\ — > 00, initially (starting from 77 = —00) grows, 
accelerating upwards, the unregularised solution diverging as r\ — > 0_. In the regularised case, 
however, in order to regain the fine-tuned behaviour as r\ — > +00, the solution has to turn around 
all the way and bend down (decelerate). This requires that there is a region around 77 = in 
which the regularised profile Q 2 {t) changes sign and becomes positive. A simple (and rather 
crude) regularisation accomplishing this can be obtained by patching the original profile outside 
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a given ^-interval to a constant potential for rj e [—770, +Vo\- In terms of the step-function H(rj) 
this means 

n 2 ( v ) = (1 - «)7 7 - 2 iT(|r 7 | - r? ) + A 2 if (r?o - M) (3.16) 

with A 2 a dimensionless constant. In terms of the original variable t the second term corresponds 
to a potential wall of width 2erjo and height A 2 /e 2 . A standard exercise in matching conditions 
reveals that for any value of the parameter a this potential admits a fine-tuned solution for an 
appropriate choice of the parameter A. 

On the other hand, the above positivity criterion immediately rules out all the candidates (3.9) 
and (3.10) since for a > 1 they are strictly negative for all r\. While there appears to be no 
analogous obvious (positivity or other) criterion for the range < a < 1, the range of principal 
interest in [26] (and also previously e.g. in [19]), it is nevertheless true that (3.9) possesses no 
fine-tuned solutions in that range either, as can be seen by explicitly solving the equations of 
motion in terms of hypergeometric functions (cf. also the comment in footnote 36 of [19]). Thus 
for this class of examples the claim in [28] (translated into our terminology) that generically 
fine-tuned solutions should exist for a discrete spectrum of values of a is confirmed in the rather 
un-useful way that this spectrum is empty. 

This shows that it is not completely obvious how to find or determine regularisations that admit 
the required fine-tuned solutions. In order to side-step this issue, we turn the question on its 
head by starting with a fine-tuned regularisation e € (t) of the solution e(t) = t 1 ~ a (and thus of 
the Rosen metric) and determining the corresponding regularised frequency f2 2 (t) , 

e(t)=t 1 - a ^e%t) => fi 2 (t) = -||| . (3.17) 

For example one can consider the even and odd regularisations 

eS(t)=t(t a + e a )-/» / ^ n e(t)-a(l-a)pn^5 ( A = 3 /( a -l) (3>18) 

Note that (3k > for a > 1 so that these regularised frequencies display the anticipated behaviour 
that they change sign for sufficiently small values of r\. Note also that for t — > these solutions 
behave as ef.(t) ~ l,t respectively, leading to the two alternative Rosen forms of the (evidently 
non-singular) Minkowski metric in that region (t = being a mere coordinate singularity in the 
latter case). 

Thus these regularised frequencies ^ 2 (t) admit fine-tuned solutions by construction. It is now 
easy to give more general multi-parameter families of regularised frequencies that accomplish 
this, simply by modifying the regularised solutions e%{t) by terms that are subleading as t — > 00 
or e -> 0. 



3.3 Adding the Dilaton 

Plane waves (more generally pp-waves) give rise to exact string theory backgrounds provided 
that they are e.g. supplemented by a null dilaton <f) = 4>( z+ ) satisfying the Einstein-dilaton 
equation R ++ = —2d + d + <f> where R ++ is the only non- vanishing component of the Ricci tensor 
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of the string-frame metric. 11 In the present (scale-invariant and isotropic) case this equation 
reduces to 

<j>(t) = 4a(a - l)t- 2 <j>{t) = -4a(a- 1) log \t\ (3.19) 

where we have selected the solution appropriate to asymptotically weak string coupling for a > 1, 
and correspondingly a strong string coupling singularity (the general solution also includes the 
general solution of the homogeneous equation, i.e. an afhne function of t). 

Once we regularise the plane wave profile, we need to enquire if the regularised dilaton equation 

<f> e (t) = -4n*(t) (3.20) 

admits solutions with these fine-tuned asymptotics or if linear terms necessarily arise that give 
rise to infinite string coupling in the far past or future. It is only in the former case that the 
pair (Q 2 ,(j) e ) can legitimately be considered to be a regularisation of the original metric-dilaton 
background. 12 

That this is indeed an issue can easily be seen explicitly in the toy model regularisation of 
(3.16) or for the regularisations given in (3.18). E.g. in the latter case the general solution of 
the dilaton equation is 

£ (t) = 2a(a - 1) [(1 - fa){t/e) arctan(i/e) - Iog(l + (t/e) 2 )] + C(t/e) + D . (3.21) 

For generic values of fa this can be arranged to have the appropriate asymptotics for either 
t — > +oo or t — > — oo by suitable choice of integration constant C (and with D = — 4a(a— 1) loge), 
but not for both simultaneously. This is only possible if fa- = 1 (with the choice C = 0), and 
comparison with (3.18) shows that this corresponds to the values fa = 1 — > a = 1 (which is 
irrelevant) or fa = 1 — > a = 4. Thus when one includes the dilaton in the story, the specific 
regularised profile given above works only for one specific value of the parameter a (and for 
other values other regularisations are required). 

The upshot of this is that the requirement that the regularised profile admit a fine-tuned dila- 
ton with weak coupling asymptotics in the past and future gives a further restriction on the 
regularised profile beyond those arising from regularity of the GY function (and its fine-tuning 
issue). We will return to the issue of how relevant this requirement actually is in the concluding 
section 4. 

3.4 Necessity to go beyond Single Scale Regularisations 

While the above seems encouraging, so far it does not allow us to claim that there is a well- 
defined evolution across t = (here we disagree, at least semantically, with [26]) and/or that 
expectation values of operators are now well-defined at t = 0. Indeed, the above analysis, in 
particular the requirement (3.13) that 3 lim e ^ (t) for ij < 0, t > does not even address the 
issue what happens at t = 0. 

11 In the Einstein frame, the equation is R\ + = (d+<j>) 2 /2 so all the backgrounds considered here satisfy the 
Einstein equations with a standard, postive-definite, lightcone energy density. 

12 The regularisation of the dilaton was also discussed in [26], but the issues that arise and are relevant in the 
present (strong coupling) context are in some sense the opposite of those encountered there. 
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Therefore let us for instance look at the e — > limit of the GY solution F t e . (t) (3.12) when one of 
the arguments is 0, say tj = 0, while t > 0. Since we have chosen the two solutions fk(v) 01 the 
regularised equation of motion (3.7) to be e- independent, the /fc(0) are finite and e-independent 
as well. One can then immediately read off that in the limit e — > or t — > oo the GY function 
is a linear combination of terms of the form e(t/e) a and e(t/e) 1 ~ a (the addition factor of e 
arising, as before, from the Wronskian) . As a consequence of the (single scale) structure of the 
regularisation, the divergence structure of the GY function for e — > is thus exactly that of the 
unregularised GY function (3.4) as U — > 0, 

e -> : F t € . =0 (t) ~ e^ 1 "", e 1 - a t a (3.22) 

and the e — > limit of the GY function F^. =0 (t) will necessarily be either or oo for alH > 0. 
Analogously, the single scale form of the regularisation implies that for t = the divergence 
structure of the would-be regularised (finite e) expectation values as e — > is identical to that 
of the expectation values (3.5) of the unregularised (e = 0) theory as t — > 0, 

e -> : {ip\z 2 \^)(t = 0) ~ Ae 2 - 2a + Be 2a + Ce 

(^|i? ; , c ^)(t = 0)^Ae- 2a +Be 2a - 2 +Ce- 1 . 

Likewise, the classical asymptotics (3.19) translate into the statement that the e — > behaviour 
of the regularised dilaton at t — is given by 4> e (t = 0) ~ — Aa(a — 1) log |e|. 

These facts illustrate that, even at the point particle level, i.e. disregarding field-theoretic infi- 
nite particle (or rather worlsheet string mode) production issues, the single scale regularisation 
advocated in [26] cannot be considered to give rise to a well-defined evolution from some initial 
time ti < across t = to some final time tf > 0. Nevertheless, as we will discuss in section 4, 
these divergences may appear in a somewhat different light when considered from the point of 
view of the space-time from which the plane wave arises as a Penrose limit. 

The above considerations do not rule out the possibility that there are more baroque regularisa- 
tions, not of the simple single-scale form (3.6), that do give finite e — > answers even at t = 0, 
and here we show, by way of example and as a proof of concept, that this is indeed possible. 
We will not elaborate on this, however, since this clearly introduces even more ambiguity into 
the regularisation procedure. 

One can for instance modify the simple regularisation (3.16) by adding a further smaller region 
around r\ = with a potential well of depth — k 2 , something like 

n 2 ( V ) = a(l - a)rT 2 # (M - Vo) + X 2 H( Vo - \r)\) - (« 2 + \ 2 )H( Vo /2 - \ V \) . (3.24) 

Analysing the conditions for existence of a GY function that is finite and non-zero in the e — > 
limit even when one of its arguments is zero, one finds that there are non-trivial solutions 
provided that k — re(e) depends suitably (and non-trivially) on e (which necessarily requires the 
introduction of new dimensionful parameters). 

Once one has a well-defined non-singular GY function F ti (t) , one can use it and its (well-defined 
and non-singular) dual solution to determine expectation values from (2.23), and these are then 
also non-singular as e — > 0. However, it is not guaranteed (and may also be too much to hope 
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for) that also the dilaton is then regularised, in the sense that the limit lim e ^ <f>e(t — 0) exists 
(and this is not true e.g. for the above profile (3.24) or simple modifications thereof that we have 
analysed) . 

4 Discussion 

Having settled, in section 2, the issue what are the right variables to use (namely Brinkmann 
variables), in section 3 we analysed various aspects of the (matrix) string dynamics in the 
singular scale-invariant plane wave backgrounds that arise naturally in this context (e.g. from 
Penrose limit [5, 6] or symmetry [13] considerations). 

It has become increasingly clear in recent years that a non-perturbative matrix string description 
in the spirit of the CSV matrix big bang model [17] is not, by itself, sufficient to resolve the 
null singularity of the space-time background (cf. also the assessment and discussion in [27]). 
For instance, as argued in [24], for the plane wave matrix big bang models of [13] near the 
singularity the usual quartic interaction term for the non-Abelian matrix string coordinates 
becomes irrelevant, due to the characteristic and inevitable presence of tachyonic mass terms 
in these models (for this reason, formally this argument applies to all the modes of the matrix 
string, not just the quantum mechanical zero mode). As a consequence, it does not appear that 
(in the spirit of the hope expressed in the context of the original matrix big bang model of [17]) 
the extra non-geometric degrees of freedom at the singularity arising from the weakly-coupled 
non-Abelian matrix string (which are actually also far from massless in the present case) can 
lead to a resolution of the singularity in this class of models. However, it would certainly be 
desirable to gain a better understanding of the non-linear dynamics of these models by other 
means. 

Moreover, as far as other degrees of freedom are concerned, the decoupling arguments of [17, 13] 
for closed and massive open string degress of freedom appear to be pretty robust. However, there 
may be some room for doubt, as the question of orders of limits becomes somewhat delicate 
when dealing simultaneously with Seiberg-Sen, near-singularity, and regularisation parameter 
e — > limits. There may also be some subtleties with the DLCQ set-up itself for time-dependent 
systems (cf. the comments in [27]). 

In the absence of an intrinsic regularisation mechanism or extra degrees of freedom in these 
models, it becomes necessary to choose some regularisation prescription in order to define the 
theory. While at this point any particular choice of prescription is necessarily somewhat ad hoc, 
one can nevertheless hope to learn something about the physics of the singularity by determining 
what kinds of regularisations are successful and what features they share. 

A very natural first step, and one which we did not question here, is to regularise the plane wave 
profile itself, as this at the very least preserves the a'-exact nature of the string background 
required for consistency of the matrix string model. Within this setting, we investigated the 
proposal of [26] to consider single-scale regularisations, depending on a single dimensionful reg- 
ularisation parameter e, which thus additionally reflect and, as much as possible, respect, the 
characteristic scale-invariance of the unregularised background. 
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As we have shown in detail, this class of regularisations can at best be considered to give a for- 
mal asymptotic (S-matrix-like) mapping between the dynamics before and after the singularity, 
without, however, providing us with a regularisation of the dynamics (time-evolution of states 
and expectation values) through the singularity. We have also traced back the failure of this 
class of regularisations to the single-scale nature of the (attempted) regularisation. 

We further found that allowing more complicated, multiple-scale, regularisations, i.e. modifying 
the profile by terms peaked around the originally singular locus and introducing additional 
dimensionful parameters beyond the overall scale e, say, one can obtain a non-singular evolution 
through the singularity. We have not elaborated nor dwelled upon this, however, because it 
evidently opens up a lot of arbitrariness and seems to be indicating that one is missing (rather 
than learning) something about the physics of the singularity. 

To better understand the physics of what is happening at the singularity in the models that 
we have analysed, recall that the plane wave metrics appearing in the plane wave matrix string 
models are to be considered as the lowest-order terms in a (covariant Penrose- Fermi [12]) expan- 
sion of the original (and singular) space-time metric (or string background) around its (singular) 
plane wave Penrose limit. Some of the divergences that one encounters appear in a different 
light from this perspective. 

In particular, one characteristic feature of the singularities of the strong string coupling (a > 1) 
backgrounds is that particles are expelled to z = oo in finite time (and strings are infinitely 
stretched, without actually going through the singularity [2, 23]). A regularisation of this be- 
haviour seems to require something novel to happen, at least some boundary condition, at 
z = oo. However, this strongly repulsive behaviour (which, see footnote 11, is not due to the 
violation of some energy condition) is not necessarily something that is exhibited by the original 
space-time metric prior to the Penrose limit, but is actually just an artefact of the Penrose limit 
truncation of the metric. 

Indeed, the Penrose limit procedure involves first a choice of null-geodesic and then a suitable 
infinite rcscaling of the metric which has the effect of infinitely expanding an infinitesimal tubular 
neighbourhood of this null geodesic to produce the entire plane wave space-time. Thus, from 
this perspective, what happens when particles reach z = oo in the plane wave is that they reach 
the boundary of an infinitesimal neighbourhood of the geodesic in finite time, which should not 
be a cause of major concern. 13 

This suggests that this particular manifestation of the divergence can be cured, or at least 
significantly altered, by retaining higher-order terms in the Penrose-Fermi expansion. This 
would also allow one to move away from null singularities, and this in turn may be good news 
anyway, as studying null singularities in a setting in which they are rigidly null may be reason 
for concern because of their potential tendency to deform into spacelike singularities. However, 
this begs the question if or how these terms can be incorporated into the matrix string theory, 
and we have no answer to this at present. 

13 Note that this interpretation is not available for the weak string coupling singularities, corresponding to the 
parameter range < a < 1, since in this case particles and strings are infinitely squeezed and excited as they 
approach the singularity. 
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Similar statements can be made regarding the regularisation of the dilaton. We investigated the 
way that the dilaton depends on the regularisation, and found that the strongest constraints 
from dilaton regularity actually come from the elimination of diverging linear dilatons at large 
times t. This would be relevant if one thinks of this as a model for big bang and emergent space- 
time, in the spirit of [17]. However, in that case one needs to study the late-time behaviour of 
the non-Abelian theory [24, 44, 45]. On the other hand if, in the spirit of the Penrose limit, one 
considers the plane wave only as a near-singularity approximation of the full metric, then the 
large-i behaviour is irrelevant simply because it is also an artefact of the approximation. 

More generally, and in conclusion, we feel that further work on these models should focus on (a) 
alternative methods to analyse and quantify the behaviour of 1+1 dimensional Yang-Mills theory 
with vanishing coupling, (b) perhaps a careful reanalysis of the DLCQ / decoupling arguments 
of [17, 13], and (c) the inclusion of higher-order terms in the Penrose- Fermi expansion in the 
matrix string models. 
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A Brinkmann and Rosen Dynamics 

A.l Classical Brinkmann and Rosen Dynamics 

The metric of a gravitational plane wave is 

ds 2 = 2dz+dz- + A ab (z+)z a z b (dz + ) 2 + S ab dz a dz b (A.l) 

in Brinkmann coordinates {z A } = {z + , z~ , z a }, and 

ds 2 = Idudv + g ij (u)dx i dx j (A. 2) 

in Rosen coordinates {x 1 } = {u,v,x 1 }. The coordinate transformation between the Rosen and 
Brinkmann forms of a plane wave metric has the form 

{u,v,x l ) = {z+,z- + \e al e\z a z\e l a z a ) , (A.3) 

where e l a — e l a (u) is a vielbein for gij(u) (satisfying the symmetry condition M ab = e a ie l b — 
ebie l a ), and the relation between gij(u) and A ab (z + ) can be compactly written as [36, 20] 

A ab (z + ) = e al {z+)e\(z+) . (A.4) 
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Following [36], it is straightforward to shows that e := dct(e") satisfies 

e/e = Tr A+ ((TrM) 2 - Tr(M 2 )) < Tr.4 = -R uu , (A.5) 

with R uu = R ++ = —5 ab A a b the only non-vanishing component of the Ricci tensor of a plane 
wave metric. In particular, if R uu > 0, then e(u ) = for some finite value of uq, and the Rosen 
coordinate system breaks down there. 

In the lightcone gauge z + — t (resp. u = t), the geodesic Lagrangian for a particle in Brinkmann 
(resp. Rosen) coordinates is 

L bc (z) = \(8 ab z a z b + A ab (t)z a z b ) , L rc (x) = y ij (t)x i x j . (A.6) 

These two Lagrangians are equivalent in the sense that under the time-dependent coordinate 
transformation x % = e l a (t)z a (the essential part of the transformation between Brinkmann and 
Rosen coordinates) they transform into each other up to a total time-derivative. 

In the body of the paper we will focus on isotropic plane waves, characterised by gij (t) = e(t) 2 5ij . 
In that case, the Hamiltonian descriptions (2.4) of the systems described by the Lagrangians 
(2.2) are related by the linear time-dependent canonical transformation 

{x,Px) = (z/e(t),ep z - ez) {z,p z ) = (e(t)x,p x /e(t) + xe(t)) (A.7) 

that can be described by a generating function of the 2 nd kind, 

{;:Sr;f. (a ' 8 > 

Given a specific Rosen metric, specified by a solution e(t) to the harmonic oscillator equation, 
one can construct its linearly independent Wronskian partner e(t), 

/"* dt' 

e(t) = e(t) / — 7^ = e(t)T(t) : W(e, e) = e~e-e~e = l , (A.9) 

Jto e \ l ) 

and consider the corresponding "dual" Rosen system, arising from the dual Rosen metric with 
9ij(t) = e(t) 2 5ij. In the classical theory, the geometric significance of the function T(t) intro- 
duced above is that of conformal time for the Rosen metric (isotropic plane waves are conformally 
flat) 

ds 2 = 2dtdx~ +e(t) 2 dx 2 = e(t) 2 (2dTdx- + dx 2 ) . (A.10) 
Its related significance in the quantum Rosen theory is explained in section 2.2. 



A. 2 Quantum Brinkmann and Rosen Dynamics 



Since classically the Brinkmann and Rosen systems are related by linear (albeit time-dependent) 
canonical point transformations, thus with a generating function that is at most quadratic in 
the canonical variables, it folllows on general grounds that these canonical transformations can 
be unambiguously implemented by a unitary (isometric) transformation in the quantum theory. 
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In the case at hand, this transformation takes the form 

Mz) = m) (z) = -^ e -^) lZ %' t {z/e{t)) 

V e M (A.ll) 

tf t (x) = (*7"Vt) (x) = y/^e-^^Me^x) . 

Here the change of argument of the wave function is just (the crucial part of) the coordinate 
transformation between Brinkmann and Rosen coordinates, the phase factor exp —ik(z, t) is due 
to the fact that the two Lagrangians differ by a total time-derivative, and the prefactor \J e(t) 
reflects the fact that Brinkmann (resp. Rosen) states are to be normalised with respect to the 
measure dz (resp. dx = dz/e(t)). 

This gives an isometry U : V. rc — > %6c of the Rosen and Brinkmann spaces of states (with 
measures dx and dz respectively) which implements the canonical transformation (A. 7) in the 
form 

U~ 1 zU = ex U^PzU = ex+— (A.12) 

e 

and intertwines the action of the Schrodinger operators id t — H , 

(idt - Hbc) M*) = -7U e - H rc W t {x) . (A.13) 



/e(t) 

To illustrate the mapping (A.ll) between Rosen and Brinkmann states, we remark that the 
Brinkmann propagator (WKB) state (2.19) remains a solution of the Schrodinger equation if 
one replaces (F,—G) by any pair of solutions (e(i),e(i) = e(t)T{t)) with Wronskian +1, and 
that this more general solution 

^ tz {z) = _!_ e (i/2)[(e(t)/e(t))z 2 -T(t)zf-2zz i /e(t)} _ (Al4) 

is precisely the one-parameter family of states one obtains by mapping the Rosen momentum 
modes ip' t k (x) (2.11) with k — —Zi to Brinkmann coordinates using (A.ll). Conversely, the 
Rosen-image of the propagator state (2.19) under U^ 1 is a momentum state of that Rosen 
metric whose metric is given by (the square of) the GY function (2.14), e(t) = F ti (t). 

In order to verify (A.13) it is useful to note that it follows from (A.12) that for any (reasonable, 
polynomial say) Brinkmann operator Obc — 0{z,p z ) one has 

U- x d b cU = 6{U- x zU, U- x pzU) - 6{ex, ex + (l/e)p x ) = (6 bc ) rc , (A.15) 

i.e. one obtains the Brinkmann operator expressed in terms of Rosen variables. In particular, 
at the level of expectation values, and for Brinkmann and Rosen states related by ip — Utp' one 
generally has 

< lP\6 b c\^> = <lp'\{6bc)rcW > , (A.16) 

while U also implements the classical relation Hbc — H rc = —dtFi (A. 8) at the operator level, 

U 1 H bc U = {Hbc)rc = H rc - (F 2 )rc 

(A.17) 

UHbcU^ 1 = (H rc )bc = H bc + {F2)bc , 
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and at the level of expectation values, 

< i>\H bc \4> > - < V'|£rc|V>' >= - < V\{h)rcW >=-< M(h)bc\i> > , (A.18) 

where (F 2 )b c /rc are the symmetrically ordered operators 



(F 2 ) bc = -\(e/e)(p z z + zp z ) + \{e 2 - cj 2 e 2 )z 2 /e 2 
{F 2 )rc = -h(e/e)(p x x + xp x ) - \{e 2 + cj 2 e 2 )x 2 . 



(A.19) 



In order to determine these expectation values for the Rosen wave packets (2.20), we need to 
calculate the expectation values of p x x + xp x and x 2 , the latter being of independent interest 
anyway. Straightforward calculations lead to 



and 



(tp'\(p x x + xpxW) = C + 2T(t)A (A.20) 
{il>'\x 2 \iP') = A + T(t)C + T 2 (t)M , (A.21) 



where 

Ak\a'{k)\ 2 B = j Akk 2 \a{k)\ 2 C = J dk ik (a(k)* a' (k) - a' (k)* a(k)) , (A.22) 

and we note for a nontrivial wave packet the terms A and B will be non-zero while C can be 
zero (e.g. if all the a(£) are real). The expectation value of z 2 can then be deduced from (2.21) 
and (A.21), 

(ip\z 2 \ip) - e(t) 2 A + e(t) 2 M + e(t)e(t)C , (A.23) 
while from (2.22) and (A.18) one finds for the expectation value of the Brinkmann Hamiltonian 

(ip\H bc \ip) = \{e 2 + Lo 2 e 2 )A + i(g 2 + u 2 e 2 )M + l(el + uo 2 ee)C . (A.24) 

Here we have made use of the fact that the dual solution e(t) (A. 9) satisfies 

e(t) = e(t)T(t) => e{t) = e{t)T(t) + e(t)- 1 , (A.25) 

and introduced the functions E(e) = \(k 2 + oj 2 e 2 ),E(e) = \(e 2 + uj 2 e 2 ) and E(e, e) = l/2(ee + 
w 2 ee). 
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